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We investigate the structure of the potential energy surfaces of the superheavy nuclei isfFmioo, 
i5gHsio8j 166H2, Jit 114, and J72I2O within the framework of selfconsistent nuclear models, i.e. the 
Skyrme-Hartree-Fock approach and the relativistic mean-field model. We compare results obtained 
with one representative parametrisation of each model which is successful in describing superheavy 
nuclei. We find systematic changes as compared to the potential energy surfaces of heavy nuclei in 
the uranium region: there is no sufficiently stable fission isomer any more, the importance of triaxial 
configurations to lower the first barrier fades away, and asymmetric fission paths compete down to 
rather small deformation. Comparing the two models, it turns out that the relativistic mean-field 
model gives generally smaller fission barriers. 

PACS numbers: 21.30.Fe 21.60.Jz 24.10.Jv 27.90. +b 



I. INTRODUCTION 

Superheavy nuclei are by definition those nuclei with 
charge numbers beyond the heaviest long-living nuclei 
that have a negligible liquid-drop fission barrier, i.e. they 
are only stabilized by shell effects [[00]. The stabilizing 
effect of the shell structure has been demonstrated in re- 
cent experiments at GSI [p|,p| and Dubna jj|, where an 
island of increased stability in the vicinity of the pre- 
dicted doubly magic deformed nucleus 2 g!]Hs || || has 
been reached. 

The full potential energy surface (PES) of superheavy 
nuclei is of interest as it allows to estimate the stabil- 
ity against spontaneous fission and to predict the opti- 
mal fusion path for the synthesis of these nuclei. Both 
features are of great importance for planning future ex- 
periments. There are numerous papers on the structure 
of the potential-energy surfaces of superheavy nuclei in 
macroscopic-microscopic models, see, e.g., [p|-|l0||, but 
only very few investigations in selfconsistent models so 
far. A systematic study of the deformation energy of 
superheavy nuclei along the valley of (3 stability in the 
region 100 < Z < 128 and 150 < N < 218 in HFB cal- 
culations with the Gogny force D\s under restriction to 
axially and reflection symmetric shapes was presented 
in O]. The full potential energy surface in the (3-j 
plane of a few selected nuclei as resulting from Skyrme- 
Hartree-Fock calculations in a triaxial representation is 
discussed in [fL2| . This investigation stresses the impor- 
tance of non-axial shapes, that lower the fission barrier 
of some superheavy nuclei to half its value assuming ax- 
ial symmetry. There is still no selfconsistent calculation 
of the deformation energy of superheavy nuclei allowing 
for reflection-asymmetric shapes. 

In a series of papers we have demonstrated the uncer- 



tainties in the extrapolation of the shell structure to the 
region of superheavy nuclei within selfconsistent mod- 
els p~3|,|l4]]. The reasons for the different behavior of 
parametrisations that work comparably well for conven- 
tional stable nuclei when extrapolated to large mass num- 
bers can be traced to differences in the effective mass 
and the isospin-dependence of the spin-orbit interaction 
fll5| . It is the aim of this paper to investigate the impor- 
tant degrees of freedom of the potential energy surface of 
superheavy nuclei for the example of a few selected nu- 



elides, i.e. 2 i«Fm 100 , ftt^ws, fi§112, ^114, and f? 2 120 



within the framework of selfconsistent nuclear structure 
models, namely the relativistic mean-field model (RMF, 
for reviews see M,0) and the nonrelativistic Skyrme- 
Hartree-Fock (SHFj approach (for an early review see 
p8|), in both cases including also reflection-asymmetric 
shapes. 



II. THE FRAMEWORK 

The comparison of the calculated binding energies of 
the heaviest known even-even nuclei with the experimen- 
has shown that the Skyrme parametri- 



tal values |13 



U. 



sation SkI4 and the relativistic force PL-40 are to be 
among the preferred parametrisations for the extrapola- 
tion to superheavy nuclei. The nonrelativistic force SkI4 
is a variant of the Skyrme parametrisation where the 
spin-orbit force is complemented by an explicit isovector 
degree-of-freedom |19|j . The energy fu nctio nal and the 



parameters are presented in Appendix A 1 . The modi- 



fied spin-orbit force has a strong effect on the spectral 
distribution in heavy nuclei and produces a big improve- 
ment concerning the binding energy of superheavy nu- 
clei @Q. The RMF parametrisation PL-40 @ aims 
at a best fit to nuclear ground-state properties with a 



stabilized form of the scalar nonlinear selfcoupling, see 
Appendix A_2 for details. It shares most properties with 
the widely used standard nonlinear force NL-Z |21|. 

Both models are implemented in a common framework 
sharing all the model-independent routines. The nu- 
merical procedure represents the coupled SHF and RMF 
equations on a grid in coordinate space using a Fourier 
definition of the derivatives and solves them with the 
damped gradient iteration method 122]. An axial rep- 
resentation allowing for reflection-asymmetric shapes is 
employed in most of the calculations, while a triaxial de- 
formed representation is used to investigate the influence 
of non-axial configurations on the first barrier. 

In both SHF and RMF the pairing correlations are 
treated in the BCS scheme using a delta pairi ng force 
p3[ Vpair = V p f n 5(ri — r 2 ), see Appendix A3 for de- 
tails. This pairing force has the technical advantage 
that the strengths V p / n are universal numbers which hold 
throughout the chart of nuclei, different from the widely 
used seniority model, where the strengths need to be 
parametrised with A dependence, and therefore in the 
description of a fission process would have to be inter- 
polated between the values for the initial nucleus and 
averaged values for the fission fragments. 

Furthermore, a center-of-mass correction is performed 



by subtracting a posteriori E CA 
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as done in the original fit of the parametrisations. 
This treatment of the center-of-mass correction is a fair 
approximation, its uncertainty for the heavy systems dis- 
cussed here is smaller than 0.2 MeV |25|| . The center- 
of-mass correction, however, has to be complemented by 
corrections for spurious rotational and vibrational modes 
as well. Their proper implementation is a very demand- 
ing task, as it requires the appropriate cranking masses. 
As done in most other mean-field calculations, we omit 
this detail. In the barrier heights, which we will discuss 
here, only the variation of these corrections with defor- 
mation enters. An estimate for these effects can be taken 
from a two-center shell model calculation of actinide nu- 
clei [£6],g7j: the amplitude of the corrections increases 
with increasing deformation, lowering the first barrier by 
approximately 0.5 MeV and the second barrier by 2 MeV. 
There is an uncertainty due to the numerical solution of 
the equations of motion which is of the order of 0.1 MeV 
even for large deformations, thus negligible in our cal- 
culations. The prescription of pairing adds another un- 
certainty to the calculated binding energies. We use the 
same pairing scheme and force for all calculations with an 
optimized strength for each mean-field parametrisation. 
The use of a local pairing force improves the description 
of pairing correlations within the BCS scheme compared 
to a constant force or constant gap approach Pq| , and 
removes some problems concerning the coupling of con- 
tinuum states to the nucleus. From possible variation 
of pairing recipes, we assume an uncertainty of the total 
binding energy of approximately 1 MeV |M§ ■ 

In the following, we will present deformation energy 
curves calculated with a quadrupole constraint (for nu- 



merical details see pOfl ), In a constrained selfconsistent 
calculation all unconstrained multipole deformations (of 
protons and neutrons separately) are left free to adjust 
themselves to a minimum energy configuration within the 
chosen symmetry. Thus the selfconsistent description of 
the potential-energy surface takes many more degrees of 
freedom into account than the 3-5 shape parameters that 
can be handled within macroscopic-microscopic calcula- 
tions. The macroscopic-microscopic models have the ad- 
ditional technical disadvantage that, for the description 
of a fission process, several nucleon-number dependent 
terms in both the parametrisation of the macroscopic and 
the microscopic model have to be interpolated between 
the values for the compound system and the fragments 
(see, e.g., ||), leading to an uncertainty of the binding 
energy in the intermediate region. It is to be noted, 
however, that there remains some open end concerning 
shapes also in the selfconsistent models as there might 
exist several local minima which are separated by a po- 
tential barrier. The numerical procedure solving the con- 
strained mean-field equations converges usually to the 
next local minimum, depending on the initial state. And 
it requires experience as well as patient searches to make 
sure that one has explored all local minima in a given 
region. 

The deformation energy curves presented in the fol- 
lowing arc shown versus the dimcnsionless multipole mo- 
ments of the mass density which are defined as 



Pl = 



47T 
3A? r 



(/ Y m ) with r Q = 1.2 A 1 / 3 fm. 



Note that these (3g are computed as expectation values 
from the actual mass distribution of the nucleus and need 
to be distinguished from the generating deformation pa- 
rameters which are used in the multipole expansion of 
the nuclear shape in macroscopic models |3l|] . Besides 
the description in terms of fit, we will indicate the various 
shapes along the paths in all figures, by the mass density 
contours at po = 0.07 fm~ . Furthermore, when looking 
at potential energy surfaces, one should keep in mind that 
these are only the first indicators of the fission properties. 
A more detailed dynamical description requires also the 
knowledge of the collective masses along the path. This 
is, however, a very ambitious task which goes beyond 
the aim of this contribution. We intend here mainly a 
qualitative discussion of the potential landscape. 



III. RESULTS AND DISCUSSION 

Figure [0 shows results of a Skyrme-Hartree-Fock cal- 
culation with SkI4 for 258 Fm, a nucleus that is located 
at the lower end of the region of superheavy nuclei. The 
strong shell effect in the prolate ground state lowers the 
binding energy of this nucleus by 19.3 MeV or 1% com- 
pared to a spherical shape, demonstrating the impor- 
tance of considering deformations for the calculation of 
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FIG. 1. Valleys in the PES of 258 Fm for SkI4 from cal- 
culations in axial symmetry with ( "refl. sym." ) and without 
("refl. asym.") reflection symmetry. In the vicinity of the 
first barrier also the result from a non-axial calculation ( "tri- 
axial" ) is shown. To give an impression of the nuclear shapes 
along the path, mass density contours at po = 0.07 fm _ are 
drawn near the corresponding curves. 



the ground-state binding energies in this region of the 
chart of nuclei. The first barrier is lowered from 11.8 MeV 
to 7.7 MeV when allowing for triaxial configurations. But 
the preference for triaxial shapes at the top of the bar- 
rier disappears when going to both larger and smaller 
deformations. It is interesting to note that the axial so- 
lutions are not continuously connected from ground state 
through first minimum when using a constraint on the 
quadrupole moment, but develop in two branches dis- 
tinguished by their hexadecapole moment. The ground- 
state branch has a diamond like shape with a /?4 much 
larger than the branch coming from outside. The contin- 
uous connection is established by the intermediate triax- 
ial shapes. 

The PES of 258 Fm shows some significant deviations 
from the familiar double humped fission barrier of the 
somewhat lighter nuclei in the plutonium region |p2| . 
There is no superdeformed minimum in the PES that 
can be associated with a fission isomer because the sec- 
ond barrier vanishes in case of symmetric breakup. This 
is due to the strong shell effect of the closed spherical 
Z = 50 shell in the two fragments, that reaches far inside 
to deformations as small as 02 = 1-0, the usual location 
of the fission isomer. This is reflected in the evolution of 
shapes along the symmetric path, that look like two inter- 
secting spheres. At large deformations around /?2 ~ 1-5 
a valley with finite mass asymmetry appears, which is 
separated from the symmetric valley by a small poten- 
tial barrier, but 5 MeV higher in energy. The occurrence 
of competing but well separated valleys and the conse- 
quences for fission or fusion complies with the results 
from macroscopic models, for a discussion see e.g. P,J33|. 
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FIG. 2. Valleys in the PES of 264 Hs for SkI4, drawn in 
the same manner as in Fig. nl 

Figure g shows the valleys in the PES of 264 Hs, at 
present the heaviest known even-even nucleus pij . Al- 
though the fragments from a symmetric breakup of 264 Hs 
are far from any shell closure, this channel of the PES 
keeps the characteristic structure of the PES of 258 Fm 
like the absence of a fission isomer and the vanishing sec- 
ond barrier. The fission path will follow the reflection 
symmetric solution, that gives a much narrower barrier 
than the asymmetric solution. Although the first barrier 
has similar width and height as the first barrier of typi- 
cal actinide nuclei, the absence of the second barrier will 
lower the lifetime against spontaneous fission dramati- 
cally. Like in the actinide region, the first barrier is a 
bit lowered if one allows for triaxial configurations. The 
reflection asymmetric solution does not lower the over- 
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FIG. 3. Valleys in the PES of 264 Hs for PL-40, drawn in 
the same manner as in Fig. hi 
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FIG. 4. Valleys in the PES of ?£jll2 for SkI4, drawn in 
the same manner as in Fig. hi 



all barrier, but it coexists far inside the barrier. This is 
a new feature occurring in the PES of many superheavy 
nuclei and was already found in macroscopic-microscopic 
calculations in the two center shell model |35|. The 
asymmetric path connects the asymptotically separated 
combination 210 Po + 54 Cr with the ground state, and cor- 
responds to the fusion path. This combination of projec- 
tile and target differs only slightly from the experimen- 
tally successful choice 207 Pb( 58 Fe, n) 264 Hs @. It reflects 
the strong shell effect of an asymmetric breakup with a 
heavy fragment in the region of doubly magic 208 Pb. It 
is interesting to compare that with the case of actinidc 
nuclei: these have also a well developed asymmetric path 
which is, however, confined to large deformations and 
reaches only down to the outer barrier | p0| . 

A calculation of the PES of 264 Hs in the RMF using 
PL-40 gives qualitatively the same results, see Fig. pi but 
there are some differences in details. The barrier is a bit 
smaller in height and width than for SkI4. The lowering 
of the barrier for PL-40 is due to the smaller shell effect 
for the ground-state configuration in this parametrisa- 
tion. While for SkI4 this nucleus has a deformed proton 
magic number, PL-40 does not predict a shell closure for 
Z = 108 at all, see |H. The effect of non-axial configu- 
rations on the height of the barrier is of the same size as 
in SkI4. 

As an example for a nucleus located at the upper bor- 
der of the known chart of nuclei, Fig. shows the valleys 
in the PES of 2 £|112, calculated with SkI4. This nu- 
clide corresponds to the compound nucleus in the cold 
fusion reaction 208 Pb( 70 Zn, n) 277 112 which was used to 
synthesize the heaviest detected superheavy nucleus so 
far H. Although the proton number of this nucleus is 
close to the value Z = 114 for the next spherical proton 
shell closure predicted by SkI4, its neutron number is 



quite far from the next predicted spherical neutron shell 
closure N — 184 but close to the deformed shell closure 
N = 162, that drives the nucleus to a strong prolate de- 
formation with /?2 = 0.22, (3 4 = —0.09, see [[t4[ for details. 
The PES of 166H2 shares most overall features with the 
PES of 264 Hs, like the one-humped structure and the 
lowering of the first barrier due to asymmetric configura- 
tions, but there are some differences in detail. The (sym- 
metric) fission barrier is narrower and slightly smaller 
(7.3 MeV compared with 10.6 MeV) than for 264 Hs. We 
have checked the effect of triaxial shapes and found that 
they are not effective to lower the first barrier for this 
superheavy nucleus. At superdeformed shapes (3% ra 0.5 
a spurious minimum develops in the symmetric barrier, 
that in reality is a saddle point, since the potential drops 
for asymmetric deformations. Note that the barrier is 
very soft in mass asymmetry in this region. Even at 
quadrupole deformations as small as fo = 0.5 the binding 
energy is nearly constant within the range < (3% < 0.3. 
The asymmetric path shows a rich substructure. There 
is a shallow minimum at fii w 0.9, while around f3% = 1.2, 
the results show a transition between two solutions with 
slightly different hexadecapole moment corresponding to 
shapes with differently pronounced "necks" but nearly 
constant mass asymmetry. The asymmetric valley corre- 
sponds to the breakup 210 Po + 68 Ni, that is quite close to 
the projectile-target combination used for the synthesis 
of this nuclide. 

Axial and reflection symmetric calculations within the 
semi-microscopic "extended Thomas-Fermi-Strutinski 
integral method" (ETFSI) |Q that uses a Skyrme force, 
i.e. SkSC4, for the nuclear interaction as well, found su- 
perdeformed minima in the PES of this and many other 
superheavy nuclei in the region Z > 112 with /% ~ 0.45 
and a larger binding energy than the usual minima at 
small deformations. Our results indicate that these su- 
perdeformed minima vanish or will have a rather small 
fission barrier when reflection-asymmetric shapes are 
taken into account. Therefore the usual minimum in the 
PES at smaller (32 has to be considered as the ground- 
state configuration, having a still sizeable first barrier 
and thus the larger fission half-life as compared to the 
competing minimum. 

All nuclei discussed so far are located in the region of 
known superheavy nuclei. Now we want to look at possi- 
ble candidates for the spherical doubly magic superheavy 
nucleus. As shown in |il|-|lS|], the predictions for doubly 
magic nuclei differ significantly between SHF and RMF. 
The RMF predicts 2 7 2 120 to be doubly magic, while the 
extended Skyrme force SkI4 prefers 2 g 8 114, the nucleus 
that has been predicted to be the center of the island of 
superheavy nuclei for a long time [|l|,|2|]. Other Skyrme 
forces, however, do not predict any doubly magic spher- 
ical nuclei in this region at all or shift the center of the 
island of superheavy nuclei to fg2l26, for example the 
force SkP |jL3| . The PES of this nucleus, calculated with 
SkP allowing for triaxial shapes, is discussed in |12] . We 
now look at the two other candidates. 




1.0 1.5 2.0 2.5 3.0 
Pi 
FIG. 5. Valleys in the PES of ?!|ll4 for SkI4 and PL-40. 
Results from the calculations in different symmetries can be 
distinguished by the mass density contours which are drawn 
near the corresponding curves. In the vicinity of the first 
barrier for SkI4 also the result from a non-axial calculation 
is shown, that lowers the barrier. 



Figure ra shows the paths of minimum potential en- 
ergy in the PES of fl! 114 : calculated with SkI4 (solid 
line) and PL-40 (dotted line). While PL-40 shows only 
a weak neutron shell closure for this nucleus, 298 114 is 
the spherical doubly magic superheavy nucleus predicted 
from SkI4. Therefore both forces lead to a spherical 
ground state of this nucleus, but with differently pro- 
nounced shell effects. In the figure, the PES from calcu- 
lations with PL-40 is shifted with respect to SkI4 in such 
a way, that the (spurious) shallow symmetric minimum 
at fa ~ 0.6 has the same energy in both models. For de- 
formations larger than fa ~ 0.5, both forces coincide in 
their prediction for the PES: The second barrier vanishes 
if asymmetric shapes are taken into account, but at large 
deformations the symmetric path is energetically favored. 
Even the shell fluctuations that lead to steps in the sym- 
metric path are located at the same deformation for both 
forces. The significant difference between the potential 
energy surfaces occurs at small deformations fa < 0.5. 
The binding energy of the spherical configuration, mea- 
sured from the reference point is lowered by 7 MeV for 
SkI4, but raised by approximately 1.3 MeV for PL-40. 
Nevertheless the spherical configuration is the ground- 
state for both forces. It remains to be noted that around 
fa ~ 0.3 the barrier is slightly lowered for triaxial shapes, 
for SkI4 by approximately 3 MeV, while for PL-40 the 
gain in energy is only a few hundred keV. 

As a final example, we consider the nucleus 172 120 
which has a spherical ground state and is a doubly magic 
system when computed with PL-40. Figure [| shows its 
PES for PL-40 and SkI4. The PES confirms the spher- 




FIG. 6. Valleys in the PES of ? 7 !l20 for SkI4 and PL-40, 
drawn in the same manner as in Fig. B. 



ical minimum for PL-40 whereas SkI4 prefers a slightly 
prolate ground state. It is to be remarked, however, that 
the actual ground state includes some quadrupole fluctu- 
ations around the minimum. In view of the weak defor- 
mation and small barrier at zero deformation it requires 
a more elaborate calculation including correlations to de- 
cide whether the true ground state will be spherical or 
deformed. A rather unexpected result is that the first 
barrier for PL-40 is indeed much lower than that for SkI4. 
This is a general result, that is also found for actinide nu- 
clei like 240 Pu. The fission half-lives from PL-40 will thus 
be generally smaller and therefore 172 120 will be more 
stable within SkI4 than with PL-40 although the lat- 
ter predicts this as a doubly magic nucleus. Both forces 
predict a strongly competing second minimum which, 
however, cannot stabilize as a ground-state configura- 
tion (or serious isomer) because the low second barrier 
makes it extremely unstable against fission. This was al- 
ready found in the previous examples and seems to be 
a general feature of superheavy nuclei. The shell effects 
cease to be strong enough to counterweight any more the 
strong decrease from Coulomb repulsion. 



IV. CONCLUSIONS 

We have presented results of constrained selfcon- 
sistcnt calculations of superheavy nuclei within the 
Skyrme-Hartree-Fock and relativistic mean-field model. 
The global structure of the PES of superheavy nuclei 
shows some significant differences compared to the well- 
known double-humped fission barrier of heavy nuclei 
90 < Z < 100. The barrier of superheavy nuclei is only 
single-humped. For the lighter superheavy nuclei we still 
find that triaxial configurations lower the first barrier, 



i.e., the region between the ground state and (3 2 »s 1.0. 
This effect vanishes for the heavier nuclei in the region 
114 < Z < 120 discussed here, but reappears in the heav- 
ier nuclei around fl®126 Mj. The second minimum (the 
fission isomer in the actinides) looses significance for all 
superheavy nuclei because it becomes unstable against 
(asymmetric) fission. Seen from the reverse side, it turns 
out that the shell structure of the final fragments influ- 
ences the PES down to small deformations. The asym- 
metric channel with 208 Pb as one fragment thus carries 
through deep into the first barrier. This corresponds 
most probably to the optimal fusion path whereas fission 
proceeds preferably along the symmetric shapes. The 
global patterns of the paths are less model dependent 
than for the actinides. Differences are most pronounced 
in the vicinity of the ground states. They are caused 
by differences in the detailed shell structure and lead to 
dramatically different predictions for shell closures and 
fission half lives. The existence and stability of super- 
heavy nuclei is thus a most sensitive probe for the present 
mean-field models. 
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and q G {p, n}. p q , T q , and 3 q denote the local density, 
kinetic density, and spin-orbit current, that are given by 



Pq= Yl v k\i>k\ 2 , T i = Y w fcl V ^ 
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Jg = -| Y V fe[^fc V X O'V'fc - (V X & Vfe) 1 'ipk 

Densities without index denote total densities, e.g. 
p = p p + p n . The ipk are the single-particle wave- 
functions and w| the occupation probabilities calcu- 
lated taking the residual pairing interaction into ac- 
count, see Appendix A3. £kin is the kinetic energy 
£kin = [fr 2 /(2m)] Jd 3 rr, while £q is the Coulomb energy 
including the exchange term in Slater approximation. 
The center-of-mass correction reads 
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APPENDIX A: DETAILS OF THE MEAN-FIELD 
MODELS 

1. Skyrme Energy Functional 

The Skyrme forces are constructed to be effective forces 
for nuclear mean-field calculations. In this paper, we use 
the Skyrme energy functional in the following form 

£ = £kin [t] + Ssk [p, t, J] + £c [p P ] - £c.m. 



with 
£sk : 



<WI |p 2 



b'r 



?E^ 



53 „<*+2 



3 



p a Y.p\ 



+b ± pr - b[ Y Pq T q - ~^P A P + i? Y Pi A Pi 



where P c . m . is the total momentum operator in the 
center-of-mass frame. The cor rection is calculated per- 
turbatively by subtracting (Al) from the Skyrme func- 
tional after the convergence of the Hartree-Fock itera- 
tion. The parameters hi and b\ used in the above def- 
inition are chosen to give a compact formulation of the 
energy functional, the corresponding mean-field Hamil- 
tonian and residual interaction p7J . They are related to 
the more commonly used Skyrme force parameters ti and 
Xi by 

b = £ (l + §0:0)5 
b' =to(j + X ), 



61 
b[ 

l>2 



\[ti{l + \x 1 )+t 2 {\ + \x 2 )]., 
J[*i(| + ari)-t2(| + x 2 )], 



I [3*i(l 



>zi 



$2(1 + 50:2)], 



b l 2 = Wh{\+x 1 )+t 2 {\+x 2 )\, 



63 = i*s(l 



5Z3) 



&3 = 1*3(3 + Z 3 ). 



The actual parameters for the parametrisation SkI4 used 
in this paper are 



to = -1855.827 MeVfm 3 , 

h = 473.829 MeVfm 5 , 

t 2 = 1006.855 MeVfm 5 , 

t 3 = 9703.607 MeVfm 3+Q , 

64 = 183.097 MeVfm 5 , 



x a = 0.405082, 

xi = -2.889148, 

x 2 = -1.325150, 

x 3 = 1.145203, 



Va 



-180.351 MeVfm 5 



with a = 0.25. For the nuclcon mass we use the value 



that gives h /(2m p ) = h /(2m n ) 
the constant entering £km- 



20.7525 MeVfni for 



2. Relativistic Mean-Field Model 

For the sake of a covariant notation, it is better to 
provide the basic functional in the relativistic mean-field 
model as an effective Lagrangian density, which for this 
study is defined as 



C 



RMF 



£ 



M 



-NM 



where 






C 



NM 



-c/ CT $ CT p s - g u $u,n(f - gp&p* 



±Ato 2 ( 6$> 2 log 



£$ 2 + (* ff - $ ) 2 



<5$ 2 + $ 2 

2<I>„ rf* 2 fc^ , ! , 



W n 



«5$2 + $ 2 y 2 

and £n is the free Dirac Lagrangian for the nucleons 
with nucleon mass ton = 938.9 MeV, equally for pro- 
tons and neutrons. The model includes couplings of the 
scalar-isoscalar ($0-), vector-isoscalar ($ w ^), vector- 
isovector (<fr P:AI ), and electromagnetic (A^) field to the 
corresponding scalar-isoscalar (p s ), vector-isoscalar (p^) 
and vector-isovector (/9 M ) densities of the nucleons. W non i 
is the stabilized self-interaction of the scalar-isoscalar 
field, behaving like the standard ansatz for the nonlinear- 
ity at typical nuclear scalar densities, but with an overall 
positive-definite curvature to avoid instabilities at high 
scalar densities [|l7|j20| . The actual parameters of the 
parametrisation PL-40 are 



du 


= 12.8861, 


m u = 


= 780.0 MeV, 


9 P 


= 4.81014, 


TO p = 


= 763.0 MeV, 


<h 


= 10.0514, 






m%, 


= 4.0 fm~ 2 , 


Am 2 = 


= 3.70015 fm" 2 , 


$0 


= -0.1 11914 fm -1 , 


<5$ = 


= 0.269688 fm" 1 



(we follow the usual convention h = c = 1 such that 
197.3 MeV = lfm -1 ). For the residual pairing interac- 
tion and the center-of-mass correction the same nonrel- 
ativistic approximations are used as in the SHF model. 



3. Pairing Energy Functional 

Pairing is treated in the BCS approximation using a 
delta pairing force J23J, leading to the pairing energy 
functional 



£» 



L?W 



d 3 r X ; 



(A2) 



where \q = - 2 Lfceo,>o fkUkVk \^k\ 2 is the pairing den- 
sity including state-dependent cutoff factors fk to re- 
strict the pairing interaction to the vicinity of the Fermi 
surface [E8[ . w 2 is the occupation probability of the cor- 
responding single-particle state and u 2 = 1 — v\. The 
strengths V p for protons and V n for neutrons depend on 
the actual mean-field parametrisation. They are opti- 
mized by fitting for each parametrisation separately the 
pairing gaps in isotopic and isotonic chains of semi-magic 
nuclei throughout the chart of nuclei. The actual values 
are 



V n 



-310MeVfm d 



V n 



-323MeVfm J 



in case of SkI4 and 



V p = -348MeVfm d 



K = -346McVfm J 



for PL-40. The pairing-active space Vl q is chosen to in- 
clude one additional oscillator shell of states above the 
Fermi energy with a smooth Fermi cutoff weight, for de- 
tails see 
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